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Abstract 

Soft- wall models in AdS/QCD generally have dilaton and scalar fields that vary 
with the fifth-dimension coordinate. These fields can be parameterized to yield hadron 
mass spectra with linear radial trajectories and to incorporate spontaneous breaking 
of chiral symmetry. We show how to construct scalar potentials which lead to such 
solutions. 



1 Introduction 



The correspondence between certain gauge theories in four dimensions and gravitational 
theories in higher dimensions has proven to be a very fruitful one. It allows for the calculation 
of physical observables in a strongly coupled gauge theory, where perturbation theory is 
inapplicable, by performing a corresponding calculation in a classical gravitational theory. 
This general gauge/gravity duality was originally inspired by the anti-de Sitter / conformal 
field theory (AdS/CFT) correspondence, and as such the metric usually takes the form of 
anti-de Sitter space in five dimensions [U [21 El II]- QCD, the gauge theory of the strong 
interactions, is strongly coupled at low energies. Heretofore, if one wanted to understand 
the structure of hadrons, then the only reliable means was via numerical computations in 
lattice gauge theory. Aside from lattice gauge theory, one must resort to models. The 
discovery of a gravitational dual to QCD would be a fantastic achievement. Unfortunately, 
the gravitational dual is not known, nor is it known whether one exists even in principle. 
Deriving such a duality from a fundamental theory, such as string theory, is referred to as 
the top down approach. The bottom up approach assumes that such a dual exists, and it 
models QCD by an effective five-dimensional gravity theory. One would like to incorporate 
the essential features of QCD into such a model. Such effective five-dimensional models are 
generically referred to as AdS/QCD, and allow for the computation of physical quantities 
in QCD such as mass spectra [51 © El El EH EI], form factors pU Q21 QU US [16], and 
thermodynamic observables [13 HH1 1201 Ell 1221 1231 

The first model to be constructed via the bottom up approach is referred to as the 
hard- wall model [HI EHJ [26J ; it simply places an infrared cutoff on the fifth-dimensional coor- 
dinate. This cutoff breaks the conformal symmetry by hand, and allows for the introduction 
of the QCD scale. The soft-wall model improves upon the hard wall, both physically and 
mathematically, by using a dilaton field to suppress the infrared contributions in a softer and 
continuous fashion [7J . This model is a mathematical improvement on the hard wall, because 
the geometry is everywhere continuous and thus avoids discontinuities and/or singularities 
which may exist in a hard-wall setup. The soft wall model is a physical improvement be- 
cause the dilaton field is parameterized as a function of the fifth dimensional coordinate in 
order to provide linear radial trajectories for the meson masses. There have been numerous 
improvements and variations on these models, some of which attempt to incorporate con- 
finement [27] and chiral symmetry breaking [28]. In the latter case, a scalar field which is 
dual to the condensation of the quark bilinear operator in QCD is included. As such, the 
behavior of this scalar field in the extra-dimensional geometry corresponds to the properties 
of chiral symmetry in the gauge theory. In an AdS/QCD model, this scalar field can be 
parameterized to incorporate both spontaneous and explicit breaking of chiral symmetry. 

As mentioned above, in bottom up models the dilaton and scalar fields are parameterized 
to reproduce certain important features of QCD, namely, confinement and spontaneous and 
explicit chiral symmetry breaking. In almost all cases, these background fields are imposed 
by hand, and are not derived as the solution to any equations of motion. There are (at least) 
two reasons that a well-defined action is desirable within AdS/QCD models. First, one needs 
a proper set of background equations in order to add perturbations to the geometry which 
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can give access to transport coefficients in the dual field theory. Second, a well-defined action 
may provide insight as to how such a model could arise from a top down approach. In order 
to meet these goals, one must know the potential which gives rise to such backgrounds. Some 
recent work on this subject was done in [29] where the authors explore the correspondence 
between the potential, and the running coupling in the dual field theory. 

A notable case where the potential was determined is the dynamical soft-wall model of 
[3D] . There, an auxiliary scalar field is introduced in addition to the dilaton, and a scalar 
potential is derived which has the soft-wall model as a solution to the equations of motion. 
Throughout the aforementioned work, some attempt is made to identify this auxiliary scalar 
field with a closed string tachyon field in string theory. Instead of this identification, one 
can promote the scalar field to a matrix valued field and identify it as dual to the quark 
bilinear operator qq which is responsible for chiral symmetry breaking. However, [30J and 
the majority of other AdS/QCD models, such as those in [7J [12j [31], have the (light) quark 
condensate proportional to the (light) quark mass, S oc m q . This is in contrast to QCD, 
where the condensate has a nonvanishing limiting value as m q — > 0. An exception is the 
phenomenological model of [2E] which parameterizes the fields in such a way as to allow m q 
and S to be varied independently. 

In this paper we attempt to bridge the gap between the dynamical model [30] and phe- 
nomenological models of chiral symmetry breaking, such as |28j, by demonstrating how to 
compute a potential given a phenomenological parameterization of the dilaton and/or scalar 
fields. This methodology is an improvement upon [30] in two respects. First, it is more 
general in that the potential can be found for many reasonable parameterizations of the 
dilaton and scalar fields. Second, even in the case considered in [30J, our resulting potential 
is both simpler and less constrained, as discussed in Sec. [4| 

Our paper is organized as follows. In Sec. [2] we identify the basic ingredients of the 
AdS/QCD model and elements of the AdS/CFT dictionary which we will use. In Sec. [3] we 
give the relevant Einstein equations, and equations for the background fields which generate 
the desired metric. In Sec. [4] we find the scalar potential when the fields have power-law 
profiles. In Sec. [5] we show how to find the potentials for models wherein the fields are 
allowed to have more complicated profiles. In Sec. [6] we provide concrete examples of some 
parameterizations which lead to analytic potentials. In Sec. [7] we discuss what happens if 
the potential is restricted to a form which has only quadratic and quartic terms in the scalar 
field. We conclude the paper in Sec. [8j 

2 Ingredients of the Model 

We are interested in dynamically generating phenomenological AdS/QCD backgrounds. Fol- 
lowing [3n] , we assume that the matter which supports the metric is a set of two scalar fields: 
(ft and x- These scalar fields interact through a scalar potential V((f>, x); it will be our goal 
to determine suitable potentials which lead to phenomenologically desirable backgrounds. 
We would like to have a background which has the following characteristics: 
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1. The metric in the string frame should be exactly five-dimensional anti-de Sitter space. 



ds 2 
Lit o q j 



string 



—dt 2 + dxidx 1 + dz 2 



(1) 



Here L denotes the AdS curvature radius, the index i runs over the three spatial 
dimensions, and z denotes the extra (fifth) dimension. 

2. At asymptotically large z the dilaton should behave as <j>{z) ~ z 2 . 

3. The mass of the dilaton should take a value such that the dual operator in the 4D field 
theory has a physically relevant mass dimension. 

Points [T] and [2] are necessary ingredients for a soft-wall model. Such a background leads 
to linear radial trajectories in the resulting meson spectrum (m 2 ~ n), in agreement with 
the data [7] . Point [3] is related to the mass of the dilaton. The relationship between the 
field's mass and the dimension A of the dual field theory operator is given by the AdS / CFT 



We will show that we are able to keep this parameter arbitrary in many of our solutions, 
but of course we want the corresponding gauge theory operator dimension to be physically 
relevant, and the mass of the dilaton should not violate the Breitenlohner-Freedman stability 
bound [32] which states that m^L 2 > —4. 

The above points are the only necessary requirements for a soft wall type of model. If, in 
addition to this, we would like to associate the scalar field x with the operator responsible 
for chiral symmetry breaking, qq, we require the following: 

4. The mass of the scalar field should be m 2 = —3/L 2 . 

5. At small z the scalar field should behave as x( z ) = Az + Bz 3 + • • •. 

6. At asymptotically large z the scalar field should behave as x( z ) ~ z - 

Because the dimension of the operator qq is 3, we require that m 2 = —3/L 2 . If the scalar 
field x is dual to the operator responsible for chiral symmetry breaking, then in the UV 
[z — > 0) regime it should have a profile which can be expanded as in item [5] above with the 
coefficients A and B proportional to the quark mass and the chiral condensate, respectively 
(cf. [U El |28]). Often in this work we will consider the chiral limit of zero quark mass; in 
this case the coefficient A vanishes, and chiral symmetry is broken spontaneously. It was 
also argued in [221 EH E2] that the scalar field x should behave linearly at large z so as to 
realize the nonrestoration of chiral symmetry at large n, again in agreement with the linear 
radial trajectories present in the datsQ Note that to properly be dual to this operator, the 

1 It is important to note there is not a consensus in the literature regarding chiral symmetry restoration in 
highly excited hadrons. For an opposing viewpoint, see |36j . The requirement of chiral symmetry restoration 
at large n would require different field profiles at large z. We will not address such a model within the context 
of this work. 



dictionary [US El El- 



m 2 L 2 = A(A - 4) 



(2) 
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field x should be charged under the bulk chiral symmetry, meaning it should be complex 
and matrix valued. We have introduced a single real scalar field for simplicity, but we will 
discuss promotion of this field to a complex, matrix valued one in Sec. [7j 

Before proceeding, let us make some comments on the relation of our approach to that 
of standard AdS/CFT. In the usual AdS/CFT approach, a conformal field theory (e.g. 
M = 4 Supersymmetric Yang-Mills theory) is deformed by introducing operators into the 
field theory Lagrangian. From the gravity point of view, one takes an existing bulk geometry 
(e.g. AdS$) and introduces fields into this background. Near the UV boundary, these fields 
act as free, noninteracting fields on AdS$. As one moves into the interior of the AdS space, 
the geometry will no longer be AdS 5 due to the backreaction induced by the new bulk fields. 

In this work, we compute all backreaction by introducing the fields into the Lagrangian 
and then solving the equations of motion for all values of the radial coordinate, including 
the interior of the AdS space. As mentioned in point [T] above, we have chosen the metric to 
be exactly AdS$, and thus the potentials which we detail in this work lead to exactly AdS^ 
metrics after all backreaction of the fields is taken into account. 

Our choice of metric is made simply to reduce the amount of freedom which exists in the 
equations of motion; the methods which we detail here could easily be extended to metrics 
which are warped versions of AdS^. Of course, one could introduce deformations (in the usual 
AdS/CFT way) into the backgrounds which we detail below ; such deformations would, in 
general, introduce backreaction which would cause the metric to be no longer exactly AdS$ 
for all values of the radial coordinate. 

Now that we have detailed the set of requirements for our model and clarified our ap- 
proach, let us now proceed to derive the relevant background equations and attempt to find 
potentials and solutions which have the characteristics we have outlined above. We will see 
that it is quite difficult to create a simple phenomenological model which satisfies all of the 
above requirements. 

3 Background Equations 

In the Einstein frame, the action can be written in its canonical form 

Se = 16^G~ 5 1 d " x ^^ ( Re ~ \ 9 ^ 9 ^ ~ \ 9 ^X ~ V(4>, X )) ■ (3) 

Here R denotes the Ricci scalar, G 5 stands for the five-dimensional gravitational constant, 
and the subscript E denotes the Einstein frame. Note that our conventions are such that 
the fields and x are dimensionless while V(<f>, x) has dimensions of energy squared. The 
energy momentum tensor which is derived from this action is 

87rG 5 T^ = | (d^d^ + d^vX - 9^C) , (4) 

C = \d x <pd x <P+\d xX d X X + V(<P,x)- (5) 

We assume that in the string frame there is a nontrivial coupling of the dilaton to the 
Ricci scalar. As such, the string and Einstein frame metrics are related by the conformal 
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transformation 



string _ 2a<f> E 
iJuu c ilu 



(6) 

where a is a constant which depends on the coupling of the dilaton to the Ricci scalar in the 
string frame action. Because the metric is static and depends only on the extra-dimensional 
coordinate z, we make the usual assumption that the fields themselves are only functions of 
this coordinate. 

There are four nontrivial background equations. We will work with the following combi- 
nations, where G^ v denotes the Einstein tensor: 



g tt G tt -g* z G zz = 8irG 5 (g u T tt - g zz T z . 

g tt G tt + g zz G zz = 8nG 5 (g tt T tt + g zz T z , 
dV 



g zz (W{z)f + [ X '(z)f 



-V(cf>, X ) 



□c 



dV 
dx 



(7) 
(8) 
(9) 

(10) 



As usual, □ = V M V M and V M denotes the covariant derivative with respect to the background 
metric. Using the presumed form of the metric expressed by Eqs. ([I]) and (|6]), these become 



"(z) + [<P\z)f (6a 2 - l) - [x\z)f + 12G0 ^ ) 



3g2a?i(z) 



L 2 



a<j)"(z) - 3a 2 [<j)'{z)] 2 - ~a<f/(z) - 



Z 



Z 

L 2 

2 r 



<j>"(z) - 3a[cf>'(z)f 



0. 

V{cf>{z),x{z)), 
dV 

<t>=<t>(z),X=x(z) 



I? 



x\z) - 3x'(z) (a(f)'(z) + 



<90 
dV 
dx 



(11) 
(12) 

(13) 
(14) 



<t>=4>(z),x=x(z) 



These equations are not all independent. Because the potential depends on z only through 
the fields, we have 

d dV ,,. . dV , 



dz 



(15) 



This relation allows one to eliminate either (14) or (|13j). 

So far we have kept the constant a arbitrary; its value can be fixed by examining (11). 
Consider this equation in the limit of large z, where the dilaton is required to be quadratic 
in z. Assuming 

(16) 



\z l 



with A a constant, Eq. (11) becomes 



\x'(z)] 2 = 36aX + 4AV 6a 2 - 1 



(17) 
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In the limit of large z, the constant term can be dropped. The solution to the resulting 
differential equation gives x( z ) ~ z<1 a ^ large z. This is in contradiction with the desired 
behavior detailed in Sec. [2| The only way to avoid this problem is by choosing 

a = ±l/v / 6 (18) 

so that the z 2 term drops out in the above equation. With this choice, one is able to have 
the desired quadratic dilaton and linear scalar field at large z. We must choose the positive 
sign so as to keep x rea l- Thus, we will fix a = 1/ v6 for the remainder of this work. It is 
quite interesting that this value of a appears frequently in the literature and can arise quite 
naturally from noncritical string theory [TQJ, [Til ESS EH EE]. With this value of the a, the 
string frame action is 

«S stri ng = / ^v^e' 2 * (R s + 4d^&>$ - IdtfPx ~ e" 4$/ V(0, X )) • (19) 

Here we are using the subscript s to distinguish the string frame from the Einstein frame. 
The field $ is a scaled version of 0, 

* = y!*- (20) 

The authors of |30j considered this action based on string theory considerations. In fact, the 
normalization of the kinetic term for the dilaton above is exactly the same as that which 
often appears in a low energy effective action for string theory (cf. Sec. 3.7 of [39]). We have 
shown that this is the only possible dilaton-scalar action which is consistent with our desired 



behavior outlined in Sec. |2| In other words, if the kinetic term in (19) were normalized 
differently, one could not satisfy all of the requirements in Sec. [2} 



4 Power-Law Solutions 



Let us first try to construct potentials which have solutions where x is purely a power-law 
in z. As such, we make the ansatz 

X (z) = XoZ n . (21) 

Of course, only certain powers of z appear in the ingredients of our model outlined in Sec. [2j 
but it is just as easy to work with a general n at this stage. Equation (11) gives the solution 
for 0, 

, . . n\/6 



12(l + 2n) AU ~ ' (22) 

Here we have assumed a Dirichlet boundary condition on so that 0(0) = 0. With this 
choice, the Einstein frame metric is asymptotically anti-de Sitter. 



Inserting the solutions for the fields into (12) and (13) yields the system 



V{(j>{z), X {z)) 



L 2 



-12+ ( * 07 ^ (2n-7) 
2(l + 2n) v ; 



(Xonz 



4(1 + 2n) 2 



(23) 
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dV 



<t>=<t>(z),X=x(z) 



(Xonz n ) 2 
2(1 + 2n 



(2n-4) 



4(1 + 2n) 2 



(24) 



As mentioned previously, we need not consider Eq. (14) since it is not independent of the 
two listed above. 

The challenge now is to determine the potential V((p,x)- The authors of [30] 



assume 

that the potential can be derived from a superpotential and use this fact to help determine 
V. We will take a different route. By examining the structure of the above equations we 
notice two facts. First, both equations contain the same exponential factor. This leads us 
to believe that the potential can be written as 



(25) 



This form of the potential is also natural on the grounds that the exponential factor arises 
due to the transformation between the string and Einstein frames. (V is the string frame 
potential.) Second, we notice that the only powers of z which appear in the above system 
of equations are z 2n and z An . Motivated by this, and the fact that we know the power-law 
behavior of both and Xi we make the ansatz 



V{<f>, X) = L2 [CO + Cl0 + C 2 X 2 + C 3 2 + C 4 X 4 + C 5 0X 2 



(26) 



The terms proportional to 
proportional to 

into the background equations (23) and (24), match the coefficients of z on each side of the 



and x gi ye r i se to the terms containing z , while the terms 
2 , x 4 ; and 4>x 2 gi ye r i se to the z An terms. One should then insert this ansatz 



equation, and try to solve for the coefficients c$...c§. 

The resulting system of equations does have a solution. Even more remarkably, the 
solution only determines five of the six coefficients. For example, C\...Cq can be written in 
terms of C3 as follows: 



Co 
Cl 

c 4 
c 5 



-12, 
4V6, 
n(n — 4) 

n 2 (c3 — 6n(l + n)) 

24(1 + 2n) 2 
n(3n — C3) 
y/6(l + 2n) ' 



(27) 
(28) 

(29) 
(30) 

(31) 



At first sight, one may be concerned about the term linear in in (26) . But when one 
expands the Einstein frame potential, one finds 



V(<f>,x) 



-12 



L 2 



4 + C3 2 n(n-4) 2 



n (4 + c 3 + 4n - 2n 2 
V6(l + 2n)L 2 



(32) 
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In fact, the value of c\ which solves the background equations is also the value which ensures 
the expansion of the potential has no linear terms. This form of the potential is exactly what 
we would expect from the AdS/CFT correspondence. The first term is the usual cosmological 
constant term, which is followed by the mass terms of the two scalar fields. Evidently, we 
should identify 

c s = <^f^ . (33) 

Some comments are in order regarding this arbitrary mass. Astute readers will notice 
that an arbitrary mass term and a field which behaves as z 2n near the AdS boundary, 
do not agree with the standard AdS/CFT correspondence where the fields near the boundary 
behave as 

<P(z -> 0) ~ Az A (l + ...) + Bz 4 ~ A (l + ...). (34) 

This discrepancy is due to two facts. First, we have two scalar fields with particular scaling 
~ X 2 near the boundary, and second there is an interaction term proportional to 0x 2 m 



(32). Consider the background equation (13) for the field near z = where the fields 



are small. The right hand side depends on dV/d(j), and using (32), we see that the leading 



term is m00 as expected. However, the final term in (32) contributes a term proportional to 



X 2 - For the solutions we consider, is proportional to x 2 near the boundary, and thus this 



term cannot be neglected with respect to the term m^cj) in (13). This is significant because 



it means that the equation of motion (13) does not reduce to that of a free scalar field on 
AdS 5 near the boundary. The original papers on AdS/CFT (cf. |4|), assumed the fields in 
question were free (i.e. noninteracting) near the boundary, and this is the reason for the 
discrepancy of our results with the usual AdS/CFT ones. If one requires agreement with 
the usual AdS / CFT prescription, one can easily do so by requiring that the interaction term 



in (32) vanish. This leads to the requirement that 

4 + c 3 = 2n 2 - 4n, (35) 

and hence 

(m^L) 2 = 2n(2n -4), (36) 

which is in agreement with standard AdS/CFT. For a quadratic dilaton (n = 1), we have 
(m<pL) 2 = —4, the same value found in [30J. Unfortunately, this requires an operator of 
dimension A^ = 2 which does not correspond to any local gauge invariant operator in QCD. 
For the remainder of this work, we will choose to continue to keep arbitrary as much as 
possible. Because the standard AdS/CFT dictionary does not apply except for particular 
values of m^, field theory interpretations of the operator dual to will require a careful 
analysis which is beyond the scope of this work. This issue does not arise for the field x- 
Regardless of the value of m^, 

m\L 2 = n(n - 4) . (37) 

Furthermore, if we require that m^L 2 = —3 as outlined in Sec. |2j the only acceptable values 
of n are n = 1 and n = 3. (This fact was also noted in the recent work |40j.) These are 



9 



exactly the two different powers of n which appear in the low and high z regions in a desired 
phenomenological model. Putting everything together, we can write the potential as 



V(<f>, X ) 



£ 2<P/V6 , 



12 + AVQ(f) + 



+ 



(m<pLy 



1 



(2 n(n — 4) 2 

? + 2 V 



(38) 



8 - (m^L) 2 + Qn\ k(j)x 2 + - [(m^L) 2 - 8 - 12n(l + n)\ k\ A 



Here we have defined 



n 



■•-12(l + 2n)- (39) 

This potential admits a solution where the string frame metric is given by 0, and the fields 
4> and x have power-law profiles (21 ) and (22). For the case of n = 1, we have the dynamical 



soft-wall model with a quadratic dilaton and linear scalar field. The potential listed above 
is an alternative to that given in [30]. We believe that this potential is an improvement over 
the latter because it is simpler and because the mass of the dilaton can be chosen at will. 

Of course, the potential listed above is not unique. One could add terms provided they 
vanish upon application of the equations of motion. If a term AV(4>, X ) is added, the above 
solution will still be a solution of the new potential provided that 



AV(<f>(z), X {z)) 
d(AV) 
d(j) 

d(AV) 



<£=<£0),X=X(z) 



dx 



0. 
0. 

0. 



(40) 
(41) 

(42) 



4>=4>( z )>X=X(z) 



where (j)(z) and x( z ) are given in (21) and (22). As an example, a term proportional to 

(0 - k X 2 Y 

will not change the equations of motion provided I >2. 



5 Solutions for General Parameterizations 



We are now in a position to find potentials for more complicated field profiles, those which 
are not exactly power laws. As mentioned previously, at small z we would like either x ~ z 
or x ~ 2 3 depending on whether the quark mass is zero or not. At large z we require X ~ z \ 
this is necessary for a quadratic dilaton (and hence linear radial trajectories) by (11). 

By examining Eq. (38 ) for n = 1 and n = 3, one notices that the only differences between 
these two potentials are the coefficients of the x 4 term and the <px 2 term. Motivated by this, 
let us make an ansatz that the potential can be written 



V{(f>, X ) 



,2<f>/V6 

L 2 



-12 + 4^0 + 



^ m4> L) 2 -8 ' 



2 -oX 2 + /i(x)x 4 + / 2 (x)0x 2 



(43) 
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The functions fx and f 2 could, in general, also depend on 0, but we have chosen them to 
only be functions of \ for simplicity. Let us also assume that x is a monotonically increasing 
function of z such that when z is large, x is large, and when z is small, \ is small. Then, to 
have the desired scalar field profile, we simply require that 



and 



/i(x -> oo) 
h(X -> 00 ) 

/i(x^o) 
/ 2 (x -> 0) 



1 

432 
a/6 
36 V 



2 mi - 32 



( 14 - L 2 m 2 



1 

432 
V6 
36 V 



2 2 



32 



f 14 - L 2 ml 



if the quark mass is nonzero. If the quark mass is zero, we desire 

3 



/i(x^o) 
/ 2 (x^o) 



784 
28 



2 m 2 - 152 



26 - L 2 m; 



(44) 
(45) 

(46) 
(47) 

(48) 
(49) 



with fx and f 2 smooth functions of \. In other words, when the fields are large, the potential 
should take the form of (38) with n — 1, and when the fields are small, the potential should 
take the form of (38) with either n — 3 or n — 1. One cannot choose any convenient 
functions fx and / 2 ; they must be consistent with the equations of motion. To see this, for 
the moment let us assume that the function x{ z ) is known and exhibits the desired low and 
high z behavior. Let us now determine the functions fxix) an d .Mx)- 

(14) using this potential 



One must go back to the original background equations (11) 
ansatz. First, one can solve (11) for (p(z) in terms of xi z )- The solution is 



1 

71 



dy 
2 



y 



x 2 [x'^)] 2 dx . 



(50) 



where 0o and 4>x are integration constants. We would like 0(0) = so that the Einstein 
frame metric is asymptotically anti-de Sitter. In order to satisfy this boundary condition, 
both integration constants must vanish. For convenience, we now perform an integration by 
parts so that 



\/6 Go 



x lx'(x)} 2 dx 



Z JO 



x 2 {x'(x)} 2 dx 



(51) 



Next, one can take a linear combination of (12) and (13) such that both the <fi" and <fi' 2 terms 
are eliminated. This equation is 



<p'(z)z + a<p(z) 



x\A 



/ 2 (x(*)) = o. 



(52) 
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where 



a: 



Substituting in the solution for <p(z) here, one can solve for / 2 



f2(x(z))x\z) 



V6 



'I -a) 



z 







x [x'( x )\ dx + a x [x'( x )} dx 



o 



(53) 



(54) 



Finally, one has to go back to either (12) or (13) and solve for fi(x)- First, let us define 
some convenient notation, 

(55) 



I( z ) = / x [x'(%)) dx 
Jo 

Then one finds the following solution for fi{x)- 



h(x(z)) x \z 
Hz 



I : - I Z X 'I M 2 I ^ 2 - - ' 



(m^L) + 3a + 



, \x{z)Y + - X \z) + ^^fi(x(z))x\z) 
2 



-f2(x(z))x 2 (z) 



a 



2(1 - a) 



-I z) 



(56) 



a/6(1 - a) 

Now we see how the functions fx and ji are correlated with the solution x{ z )- 

Of course we are not finished, because the potential should be a function of the fields 
0, x only and should not depend explicitly on the coordinates. The expressions for fi and 
/2 above need to be rephrased so that they only depend on the field x- To achieve this, 
instead of parameterizing x{ z )i we should specify the inverse relationship z(x). As before, 
we assume that z is a monotonically increasing function of x- F° r the correct asymptotic 
behavior, we require 

z(x^°°)~X, (57) 



and 



if the quark mass is nonzero, and 



z{x -> °) ~ X 



z(x -> 0) ~ x 



1/3 



(5f 



(59) 



if the quark mass is zero. It is now a simple matter to transform the potential using the 
relation 

X\z) = ~ . (60) 



For example: 



i{x) 



z'(x) ' 

xO) 1 dz 

x(o) [z'(x)\ dx 



x Ax) 

o z'{x) 



dx 



(61) 
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This allows us to compute the potential as a function of the fields only. For convenience, 
define 



Ux) 
6(x) 



a 



1-aJo z'(x') 

l rx z(x'^ 2 



z(x) Jo z '(x') 
Then the solutions for and the potential are 



<{>(X) 



x 2 / 2 (x) 



1 



1 — a 



x 4 /i(x) 



In all, 



i (z(xY 

" 2 \z'( X ), 

1 ( z(x) ' 

2 U'(x). 



" V6 
1 — a f 



a 



a(x)-6(x) 



+ 2* + 4a 



a)[6(x) + 6(x)] 



a - l)&(x) + 6(X)] 2 + ilix) ~ 166(X) 



L 2 

\ 2 



-12 + ^0[8+(l-a) (6(X)+6(X)) 



3a 



+ 



a — 1 

4a 



a)[£i(x)+6(x)] 2 -£ 2 2 (x) + i6£i(x) 



(62) 
(63) 



(64) 
(65) 
(66) 



(67) 



This set of equations is one of the central results of this work. Given a phenomenological 
parameterization z(x), one needs to do two integrals to determine £i and £2, after which 
point the potential which gives rise to the desired solution can be determined. In addition, 
one can see how the solution for must be correlated with such a parameterization from 



(64). 



The simplicity of the potential is dependent on the simplicity of the £ functions, which 
are in turn related to the parameterization z(x)- We are free to choose z(x) at will. Unfor- 



tunately, there is no guarantee that the potential has an analytic form, as the integrals (62) 



and (63) cannot always be done in closed form. In the following section, we will give a few 



examples of parameterizations z(x) which lead to analytic potentials. 



6 Examples 



A simple parameterization z(x) which has the correct behavior when the quark mass is zero 
is 



z{x) = 7 



X+[f 2 



1/3' 



(68) 
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where f3 is a positive dimensionless constant and 7 is a positive constant with dimension 
of length. Such a parameterization allows for an analytic potential. The relevant functions 
which appear in the potential are 



Ux) 



4 



a 



27/32 (yi/3 + 2,) 



I8f3 2 (a- 1) 
1/3 



6y^-9y^-9y>-2ln(l + 3y^)} , 



y--y + -y 5/3 + —y 7/3 
y y 15 y 28 y 



27 , 
+ I2 S 



73 



arctan {\f2>y 



1/3 



where y = fix- in fact) the above parameterization can be generalized to 



z{x) = 7 



X_ 

J 2 . 



1/3 



1 + iPx) 



2/3ra 



(69) 
, (70) 

(71) 



Both of the integrals involved in the computation of £1 and £2 can be done analytically 
without the use of special functions if n is an integer (though n must be positive in order for 
the field x to have the correct asymptotic behavior). The complexity of £1 and £ 2 appears 
to increase with increasing n, so we have quoted the simplest example above with n — 1. 
Another parameterization for zero quark mass which leads to the desired behavior is 



z{x) 



IX 



1/3 



/3 2 / 3 {1 + (/3x) 1/3 [arctan ((f3 X ) 1/s ) - vr/2]} 



(72) 



This parameterization was found by examining the behavior of the integral appearing in 
£ 2 (x)- The integrand must behave as x^ 3 f° r small x> and as x 2 f° r large x- Solving the 
differential equation 



z\x) 



7 



(/?X) 4/3 + iPxf 



(73) 



z'ix) P 2 

leads to the above parameterization. With this parameterization, the solutions for £1 and £2 
are 



6(x) 



-y 2/3 + \y i/3 + 6 jy 2 + 7y 8/3 



a 

7p 2 (l - a) 

+ (jy 3 + 9y 7 ' 3 ) (arctan^ 1 / 3 ) - fj + ln(l + y" 3 ) 



Ux) 



9y 2 + 7y 8 / 3 
7f3 2 



1 + y 1 ^ 3 ( arctan(y 



l/3\ _ 



7T 



(74) 



(75) 



where again y = (3x- This can be generalized to other potentials by solving the differential 
equation 



z 2 (x) 
z'(x) 



7 d 
Wdx 



Wx) 7/3 Q + {fix) 



2/3n 



(76) 



and requiring z to have the correct asymptotic behavior. This method again leads to an 
analytic potential for n being an integer, though the simplest result is that given above with 
n= 1. 
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It is tempting to try to use functions such as the exponential and hyperbolic tangent 
to parameterize z(x). However, we have not found any applicable parameterization using 
these functions where both of the relevant integrals appearing in £1 and £2 have an analytic 
solution. 



7 Potentials Quadratic and Quartic in the Scalar Field 



As mentioned in the Introduction, in order for the field x to be dual to the operator qq it 
should be a complex, matrix valued field. So far, we have only considered x to be a real 
scalar field. In order to address this difficulty, one can simply promote the field with the 
replacement 

(77) 



:X 2 -> X f X, 



where the extra factor of 1/2 is introduced to give a canonical action for a complex scalar 
field. However, one will notice that all of the potentials which were constructed in the 
previous section contain fractional powers of Xi an d hence such a replacement is less than 
desirable. If the action is a function of x 2 om y> then such a promotion is possible. Below, 
we will discuss two methods which allow one to construct potentials which depend only on 
X 2 

The first method we use to determine a potential which is a function of x 2 only is to 
modify the ansatz (43 ) so that the functions f\ and f'2 are functions of <fi only. More generally, 
let us assume 



V(<f>, X ) 



L 2 



F G {<j>) + F 2 (0) X 2 + F 4 (0) X 4 



(7E 



with the F functions to be determined. This strategy involves parameterizing z(4>) instead of 
z(x), but the general steps involved are similar to the cases already discussed. The challenge 
is to find a parameterization with the desired behavior that leads to an analytic expression 



for x(4>)- The latter can be determined from (11) to be 



X(0) = ±6 1/4 



d 



1 1/2 



In 



(79) 



In practice, it often is easiest to start somewhere in the middle by parameterizing some 
intermediate quantity such as x' (</*)• ft turns out that a parameterization of the form 



[x'(0)] s 



1 + 



30) 



has the desired asymptotic behavior and leads to analytic solutions for both x{4>) an d z(<f>). 
In general, these solutions involve hypergeometric functions. A notable special case is if 
n — 1/2. Then the solutions are 

7 0l/6 



6 



(/V0) 1/3 -(i + /V0) 1/3 



1) 
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x(<P) 



2 3/4 



G{<f>) - G(0) + In 



8 + 3^ 



+ 9/?v^ + G(0) 



with 



G(<f>) = 3^(1 + /3 v / 0)(7 + 9/3^) 
For this parameterization define 



2 z'(0) ' 



and make use of the calculus relations 

0"(z) = 
x'W = 

x"0) = 



-A0)M0)] 3 , 

x '(<f>V(<f>) 



(82) 



(83) 



34) 



(85) 
(86) 

(87) 



Then the background equations ( 12 ) and ( 14 ) can be written 



W) + W)x 2 (</>) + W)X 4 (0) = -12 - 8R(<f>) 



2F 2 (0) X (0)+4F 4 x 3 (0) = ^X"(0)^ 2 (0) 



R 2 {<j>) + \R\<P) [ X \<P)f 



R{4>) - -R{(j>) [x'm 



9) 



3" vr/ ^' 3 

These equations can be solved to determine the F functions, and hence the potential in 
terms of the known functions xi'P) an d R{4>)- Notice that there are only two equations for 
the three F functions; we have some freedom to choose one of the F functions at will. One 
possible choice is F = —12 + Ay/6(j) — |a</> 2 , which is the same as the form in the ansatz 
(43). The resulting solutions for F 2 and F 4 are quite complicated. In addition to this fact, 
the leading term in the small field expansion of _F 4 is -5 / 6 , which is in contradiction with a 
well-defined conformal limit: V(<f) — > 0, x ~ * 0) = — 12/ L 2 . 

On the basis of algebraic simplicity, an alternative is to choose F 4 to be constant. It is 
then straightforward to solve for F and F 2 . For small values of 0, these functions have the 
expansions 



F o (0) = -12 + — x 



4 + 3x 1/6 + 0(x 1/3 ) 



(90) 



9 



.1/6 



+ 0(x 



1/3^ 



(91) 



3 

"2 ~ 2" 

with x = f3 2 <f). This solution appears to be consistent; however, it suffers from the drawback 
that the potential is quite complicated, and that the small field expansion of F contains 
fractional powers of <p starting with 7//6 . Thus, in this solution the mass of the dilaton is 
not well-defined. This solution adheres to all of the ingredients we set out in Sec. [2] except 
for point |3j 
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We now detail a second way to determine a potential which is a function of x 2 an d is 
thus a good candidate for a dynamical model of chiral symmetry breaking. As above, we 
will sacrifice point [3] in the ingredients of the model by choosing a nonstandard dilaton mass. 

One may have noticed that the general potential (67) simplifies greatly in the special 
case a = 1. This corresponds to m^L 2 = 5. If one naively uses the AdS/CFT dictionary to 
compute the dimension of the corresponding operator in this case, one finds = 5, which 
is nonrenormalizable. We are unsure whether it is simply coincidence that this choice 
simplifies the potential greatly, or whether there is some physics hidden here. Such a choice 
is certainly nonstandard, though it could be acceptable within the context of effective field 
theories. However, it should be emphasized that in light of the discussion in Sec. |4j any use 
of the standard AdS/CFT dictionary in regards to the field <fi should be approached with 
caution due to the fact that it does not reduce to a free field near the AdS boundary. 

With the choice of (m^L) 2 = 5, the potential (67) become^] 



V{<j>, X ) 



~1X~ 



j-12 + 4^60- ^0 2 + 



v/6 



<t>Kx) 



\I{X) 2 -±I{X) + \[I\X)?\, (92) 



where I(x) is defined in (61). Notice that only one of the integrals appears, and all reference 
to the function £2 has disappeared in the potential. At this point we can simply parameterize 
I(x) to fit our needs. Any parameterization of I will do, provided that for small or large \ 
the function behaves as 

(93) 



n 



Hx) -> ^X' 



where n is the desired power of the scalar field in this regime (i.e. x{ z ) ~ z n )- F° r example, 
if one chooses 

I(x) = t ( 3 + V>X)' \ (94) 

with (3 a constant, one gets the desired asymptotic behavior for zero quark mass. All factors 
of x 2 can De promoted to x^X without any problem. Another such parameterization is 



Kx) 



* \ 1 1 



1 + P 2 x 2 



(95) 



The solutions for the fields <p(z) and x( z ) could then be found numerically from (61) and 

del. 



8 Conclusion 

In this paper we addressed the problem of constructing models of AdS/QCD. Specifically, we 
showed how to construct a potential for dilaton and scalar fields that leads to an AdS/QCD 
model with many essential features of QCD. Given a suitable parameterization z(x) or z(<j)), 

2 One should take care when making this simplification due to the presence of (a — 1) in the definition of 
the function £1. 
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we show how to construct a potential V((f>, x) which will have a solution with the desired 
properties; as such, the main results of this paper are (64) and (67). Linear radial trajec- 
tories, conformal symmetry breaking, and both spontaneous and explicit chiral symmetry 
breaking can be incorporated. The desired ingredients of the model as detailed in Sec. [2] 
are numerous, and the fact that such a solution can be dynamically generated at all within 
such a simple setup is somewhat surprising. It is especially interesting that the mass of the 
dilaton can be kept arbitrary throughout much of the analysis. However, as discussed in 
Sec. |4j traditional field theory interpretations of dual operators via the standard AdS/CFT 
dictionary are questionable, except for particular values of the dilaton mass. Explicit ex- 
amples of potentials were given, at least for the case when the light quark mass is zero, 
although this was only for purposes of illustration and not a limitation in principle. We 
have also not discussed the stability of the examples we provide, as that analysis would be 
beyond the scope of this work. Before calculating physical observables within the context 
of such a solution, its stability would have to be checked. In order to incorporate chiral 
symmetry breaking, the scalar field we introduce into our action needs to be promoted to a 
matrix valued field, and this requirement complicates the analysis. However, we find that a 
nonstandard dilaton mass choice greatly simplifies the results. 

One might have expected to find simple expressions for the potential since the Lagrangian 
for QCD is so simple to write down. However, there is no reason for this to be so and, un- 
fortunately, the analytic expressions in our illustrative examples are not so simple, although 
the general structure is. It remains a challenge to find a simple expression for the potential 
that leads to the desired properties for the dilaton and scalar fields. Nevertheless, the poten- 
tials given here can be useful in the context of hadronic structure. We also believe that the 
methods we have outlined here may be useful for the determination of potentials in other 
bottom up AdS/QCD models, and at finite temperature. 
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